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Abstract. The upcoming Karlsruhe Tritium Neutrino (KATRIN) experiment will put 
unprecedented constraints on the absolute mass of the electron neutrino, . In 
this paper we investigate how this information on nii,^ will affect our constraints on 
cosmological parameters. We consider two scenarios; one where to^^. — (i.e., no detection 
by KATRIN), and one where m^^ = O.SeV. We find that the constraints on m^^ from 
KATRIN will affect estimates of some important cosmological parameters significantly. 
For example, the significance of < 1 and the inferred value of Q\ depend on the results 
from the KATRIN experiment. 
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1. Introduction 

The large amount of new cosmological data in the last decade has lead to what one 
may call the cosmological standard model. In this model the universe is close to flat, 
homogeneous and isotropic on sufficiently large scales, and today the energy density of the 
universe is dominated by dark energy (~ 74%), dark matter (~ 22%) and baryonic matter 
(~ 4%). This model is consistent with data ranging from the WMAP measurements of 
the anisotropics of the cosmic microwave background (CMB) radiation [1] to observations 
of supernovae of type la, galaxy distributions and several other observables (with a few 
exceptions, see [2]). It is often claimed that most of the data can be fltted with only six 
free parameters. This claim rests on the assumption of massless neutrinos, an assumption 
justified by the fact that adding the sum of the neutrino masses as a free parameter does 
not improve the fit substantially. 

However, from the observation of neutrino oscillations, there is a compelling body of 
evidence for non-zero neutrino masses (see [3] for a review). Oscillation experiments do not 
give us any information on the absohite mass scale of neutrinos, only on the mass differences 
between the different mass eigenstates and mixing angles. The current best upper bounds 
on the neutrino mass from particle experiments come from the Troitsk [4] and Mainz [5] 
tritium beta decay experiments that found upper bounds on rrii,^ of rrii,^ < 2.2eV (95% 
C.L.). The KATRIN experiment [6] that will start taking data in 2010, is expected to 
lower this limit on m^^ by an order of magnitude to m^^ < 0.2eV (in the case of no 
detection) after three years of running. 

Effects of neutrino masses can also be seen in cosmological observables, and the best 
upper limits on the absolute scale of the neutrino mass today come from cosmology. Both 
CMB and the large scale structures (LSS) of the galaxy distribution are probes that are 
sensitive to the neutrino mass, the observable quantity being the sum of the three neutrino 
mass eigenstates, M^, = J2i^u,- The upper bounds on from cosmology range from 

< 0.2eV [7] to < 2.0eV [1] (95% C.L.), depending on the data [8] and cosmological 
model [9] used. 

On the experimental side there is a claim of detection of the absolute scale of the 
neutrino mass from the Heidelberg- Moscow neutrinoless double beta decay experiment, 
with an effective electron neutrino mass of {m^^) — (0.2 — 0.6)eV (99.73% C.L.) [10]. 
However, these results are regarded somewhat controversial. The cosmological implications 
of this result are discussed in [1 1] . 

We know that neutrinos are massive, and since we have no current priors on the 
neutrino mass in the allowed cosmological range, one should always marginalize over M^, 
when constraining other cosmological parameters. Mi, turns out to be partially degenerate 
with several of the standard cosmological parameters, such that this marginalization over 
Mj, weakens the bounds on the other parameters in our model. Thus, any prior knowledge 
of Ml, from non-cosmological experiments will serve to break degeneracies and improve 
the constraints on other cosmological parameters. The KATRIN experiment will provide 
us with such a prior on Mi, in a range that is relevant for cosmology. In this paper we 
investigate how the results from KATRIN will affect our estimates of other cosmological 
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parameters. 

Limits on the neutrino mass when combining resuhs from KATRIN and WMAP have 
been studied in a recent paper by that H0st et al. [12]. Our emphasis in this paper is on 
how other cosmological parameters are affected when the results from KATRIN are used 
as an external prior. 

Section 2 contains a short review on the effect of massive neutrinos in cosmology. In 
section 3 wc will present the data and methods that we will use in our analysis, including the 
assumed priors from the KATRIN experiment. Then we will present our results in section 
4. A comparison of values found when introducing the KATRIN priors is presented in 
section 5. Finally we summarize and conclude in section 6. 



2. Cosmology w'lih massive neutrinos 

All our results are derived within the standard cosmological paradigm of a flat ACDM 
model, using the following free parameters: {Qft/i^, Q^, log(10^°y45), /i, n^, r, M,^}. Here fij 
denotes the energy density of energy component i (m=matter, b=baryons, A=cosmological 
constant, CDM=coId dark matter) relative to the total energy density of a spatially flat 
universe. The matter density, n„j, is the sum of all non-relativistic components, such that 
Vtm = ^CDM + fib + ^i/- The parameter h is the dimensionless Hubble parameter, defined 
by Hq = 100/ikm s~^Mpc~^, Ag denotes the amplitude of the primordial fluctuations, 
while Ug gives the tilt of the primordial power spectrum. Finally, r is the optical depth 
at reionization. For more details on the parameter deflnitions, see the description of 
the CosmoMC code [13]. The effect of massive neutrinos on cosmological observables is 
parameterized by M^,, the sum of the neutrino masses, and is related to the neutrino energy 
density by the simple relation [14] Q^.h'^ = gg • ^^^^ ^^^'^ extend the parameter space 
by including w, the equation of state parameter of dark energy, as a free parameter. We will 
assume w to be constant. This parameter may be interesting to study, as it is fundamental 
in the understanding of the nature of dark energy, and since it is known to be slightly 
correlated with [15]. It should be stressed that this analysis rests on the assumption 
of a standard thermal background of 3 weakly interacting neutrino species. Alternatives 
to this picture are studied in e.g. [18] and [19]. 

Recent reviews of the role of massive neutrinos in cosmology can be found in [14, 16]. 
In this section we will only give a brief description of the most important effects of M„ 
on relevant cosmological observables. We will throughout this work assume that the 
neutrino mass eigenstates arc degenerate, such that Miy = 3m^^^. In the mass range that 
we will operating in, it has been shown that this is a valid simpliflcation when it comes to 
cosmological observables [17]. 

Effects on the CMB from massive neutrinos manifest themselves mainly on the level of 
background evolution. In the neutrino mass ranges relevant to us, the neutrinos will still be 
relativistic at the time of matter-radiation equality, and must be regarded as a radiation 
component when it comes to the background evolution of the universe. Increasing M^, 
(and thus also fl,^), keeping flm constant, will thus postpone the time of matter-radiation 
equality. This will enhance the acoustic peaks in the CMB power spectrum and give a 
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Figure 1. CMB (left panel) and LSS (right panel) power spectra with different values of 
Ml,. Here Vl„i is held constant, and increasing Af^ has been compensated by decreasing 
D M correspondingly. 



small horizontal shift of the peaks to larger scales. This effect is shown in the left panel of 
Figure 1. To compensate for this effect, one can increase VLm and decrease Hq. It is already 
obvious that My will be correlated with both Hq and Vt^ (and thus also when we stick 
to the assumption of spatial flatness). Another effect comes from neutrino free-streaming, 
which will smoothen out gravitational wells on scales below an M^- dependent neutrino 
free streaming scale [20,21]. On scales smaller than this, the acoustic oscillations will be 
enhanced, increasing the height of the peaks in the CMB power spectrum. 

Neutrino masses affect the LSS power spectrum in an even more distinct way. Again, 
massive neutrinos will suppress structure growth on scales below a free streaming scale 
given by [14] 

knr = OMO^l^^hMpC-\ (1) 

The smaller M^, the larger scales will be affected, and the larger M^y, the more suppression 
of power on the scales affected. The effect of neutrino mass on the matter power spectrum 
can be seen in Figure 1. Again, is kept constant, and increasing M^, has been 
compensated for by decreasing Qcdm correspondingly. 



3. Data and methods 

3.1. Cosmological data 

Our analysis include both observations of CMB, LSS, SNla, information about baryonic 
acoustic oscillations (BAO) in the matter power spectrum and constraints from the cluster 
mass function from weak gravitational lensing. We have also applied priors on Hq and Qf,. 

The CMB data used in our analysis, comes from the temperature [22] and polarization 
[23] data from the 3-year data release from the WMAP team. The WMAP experiment is 
a satellite based full-sky survey of the CMB temperature anisotropics and polarization. In 
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our analysis of the WMAP data we have used the Fortran 90 hkehhood codej provided 
with the data release. 

We have used LSS data from the Sloan Digital Sky Survey (SDSS) luminous red galaxy 
(LRG) sample [24]. As SNla data we have used the sample from the Supernova Legacy 
Survey (SNLS) [25] . Other probes of the matter distribution that we have applied come 
from the measurement of the baryonic acoustic peak (BAO) in the matter power spectrum 
and the cluster mass function (CMF) from weak gravitational lensing. The BAO constraint 
comes from the SDSS-LRG sample [26], and we have adopted the fit function from [27], 



^ N-0.35 



^BAO = 0-469 [-^) (1 + 0.94/,) ± 0.017, (2) 

where 



^BAO 



Dm{z? 



^ ° , (3) 



and Dm{z) is the comoving angular diameter distance at redshift z. 

Handles on parameters governing the clustering of matter are also provided by the 
cluster mass function. The cluster mass function from weak gravitational lensing, as 
measured in [28] , gives constraints on a combination of Q^n and cTg (the root-mean-square 
mass fluctuations in spheres of radius 8/i~^Mpc). We have adopted the fit-function for 

XCMF = 10000"^ + 6726m^ + 1230^2 - 4.09m + 0.004 (4) 

where u = asiVlm/O.Sf'^^ - 0.67. 

A prior on the Hubble parameter, h = 0.72 ± 0.08 [29] comes from the Hubble Space 
Telescope (HST) key project. From big bang nucleosyntesis (BBN) we adopt a prior on 
the baryon density today, fi^/i^ = 0.022 ± 0.002 [30-32]. 

Throughout the entire work we also apply a top-hat prior on the age to of the universe: 
10 Gyr < < 20 Gyr. 



3.2. Constraints from KATRIN 

The KATRIN [6] experiment measures the energy distribution of electrons from tritium 
beta decay. The exact shape of the end of this spectrum will depend on how much of 
the energy that is bound in the outgoing electron neutrinos, and thus also be a probe 
of the electron neutrino mass. If KATRIN does not detect m,^^, they are expected to 
place an upper limit on rriu^ < 0.2eV (90% C.L.). They expect to reach an uncertainty of 
(7^2^ ^ 0.025eVl 

Here we have adopted this uncertainty for two cases, one assuming m,,^ = OeV (i.e., 
no detection by KATRIN), and one assuming m^^ = 0.3eV (giving M^, = 0.9eV). Further 
we have assumed the Gaussian distribution of ml^ around these values [6,12], and used 
this as a prior in our cosmological parameter analysis. 

I http://lambda.gsfc.nasa.gov; version v2p2 
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3.3. Method 

Employing the publicly available Markov Chain Monte Carlo code CosmoMC[13] we have 
studied our seven-parameter model for two combinations of datasets; first using only 
WMAP data, and then adding LSS, SNla data and priors from HST, BBN, BAO and 
CMF. In both cases we have compared the results from using only cosmological data, and 
from adding priors from KATRIN in the case of m,^^ = OeV and m,^^ = 0.3eV. First we will 
assume w = —1 (cosmological constant). 

Yet more freedom in the cosmological model might be added by including w as a 
free parameter, yielding a more general form of the dark energy component. We will also 
include w in our analysis, assuming it to be constant. 

4. Results 

4.1. A 7 -parameter model 

Starting out, we considered the simplest case using only the standard 7 parameter universe 
model, as explained in section 2, and WMAP data only. Then we added the assumed 
KATRIN priors for m,^^ = OeV and m^^ = 0.3eV as explained in section 3.2. The results 
are summarized in Figure 2 and Table 1. One easily sees that the different priors from 




Figure 2. Marginalized parameter distributions when using WMAP data (black solid 
lines) , compared to the resulting distributions when adding KATRIN data with m^^ — OeV 
(red dashed lines) and niij^ = 0.3eV (blue dotted lines). 
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Tjiv^ = free m,^^ = OeV rrij,^ = 0.3eV 





0.0216 ± 0.0008 


0.0220 ± 0.0007 


0.0215 ± 0.0007 


h 


0.651°:°^ 


0.70 ± +0.04 


0.63 ±0.03 


T 


0.085 ± 0.029 


0.089 ± 0.029 


0.085 ± 0.028 


ris 


0.940 ± 0.022 


0.956 ±0.016 


0.939 ±0.016 






0.71 ± 0.06 


0.60 ±0.05 




0.66 ±0.07 


0.73 ±0.04 


0.65 ±0.05 




0.34 ±0.07 


0.27 ±0.04 


0.35 ±0.05 




11.2 ±2.6 


11.3 ±2.6 


11.2 ±2.7 


(eV) 


< 1.75 (95%C.L.) 


< 0.58 (95% C.L.) 


0.871°:}° 



Table 1. Limits on cosmological parameters with different priors on the neutrino mass 
when using WMAP data only. In the left column are the results from having no priors on 
My (black solid lines in Figure 2), the middle column shows the results when using the 
assumed KATRIN prior in the case of m^^ = OeV (red dashed lines in Figure 2), and the 
rightmost column gives the results with an assumed KATRIN prior for m^^ = 0.3eV. 

KATRIN indeed affect some of the parameter constraints. 

Evidently the inferred vahie of h depend heavily on which prior we assume on m^^. A 
larger My requires a smaller h. This degeneracy is not very surprising, as both Afj, and Hq 
tend to shift the positions of the acoustic peaks in the CMB power spectrum, as mentioned 
in section 2. Note that we have assumed no prior on Hq in this case. Also erg is highly 
dependent on the priors on rrii,^ . We saw in section 2 that alters the height of the peaks 
in the CMB power spectrum, and it will therefore be strongly correlated with erg which is 
an amplitude parameter. 

Effects on Ug may be even more interesting. The significance of < 1 is important, 
as rig ^ 1 is a generic prediction of most inflation models. In the case of rriv^ = OeV, we find 
a significance of 2.7a for rig < 1, while this increases to 3.8a when applying a KATRIN 
prior in the case of rriy^ = 0.3eV. This degeneracy between and Ug stems from the fact 
that increasing tends to increase the amplitude of the acoustic peaks on scales smaller 
that the neutrino free streaming scale. This can be compensated by decreasing Ug which 
will give a larger tilt on the primordial power spectrum. 

Modifications in the distributions of ftrn a-nd are also evident. This happens since 
increasing M^, will shift the time of matter-radiation equality and thus amplify the acoustic 
peaks, and this effect can be compensated by increasing (and thus also reducing ^a)- 

Next, we added data from LSS, SNl, HST, BEN, BAG and CMF to the WMAP data. 
Doing this, the limit on from cosmology alone is in the same range as we get from 
KATRIN in the case of rriy^ = 0, thus we would not expect a large effect from adding the 
KATRIN prior in this case. If, however, we apply the = 0.3eV scenario, we would 
expect to see some effects. This is indeed the case, as we can see from Figure 3 and Table 
2. 

The effect of adding a KATRIN prior with m^^ — 0.3eV is most pronounced for the 
^m, CTs and h. The shifts in the distributions can be explained by much of the same 
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Figure 3. Marginalized parameter distributions when using data from 
WMAP+LSS+SNl+HST+BBN+BAO+CMF (black solid lines), compared to the 
resulting distributions when adding KATRIN data with niy^ — OeV (red dashed lines) 
and rriy^ — 0.3eV (blue dotted lines). 



rriy^ = free m,y^ = OeV m^^ = 0.3eV 





0.0220 ±0.0006 


0.0220 ±0.0006 


0.0221 ±0.0007 


h 


0.70 ±0.02 


0.71 ±±0.02 


0.68 ±0.02 


T 


0.086 ±0.029 


0.084 ±0.028 


0.100 ±0.029 


Us 


0.956 ±0.014 


0.955 ±0.014 


0.958 ±0.015 


0-8 


0.70 ±0.03 


0.71 ±0.03 


0.66 ± 0.04 


^^A 


0.74 ±0.02 


0.74 ±0.02 


0.71 ±0.03 




0.26 ±0.02 


0.26 ±0.02 


0.29 ±0.03 




10.9 ±2.6 


10.8 ±2.6 


12.3 ±2.5 


M, (eV) 


< 0.55 (95%C.L.) 


< 0.47 (95% C.L.) 


0.58l°i? 



Table 2. Limits on cosmological parameters with different priors on the neutrino mass 
using the fuU range of data sets WMAP+LSS+SNla+HST+BBN+BAO+CMF. In the 
left column are the results from having no priors on Af^ (black solid lines in Figure 2), 
the middle column shows the results when using the assumed KATRIN prior in the case 
of rriu^ = OeV (red dashed lines in Figure 2), and the rightmost column gives the results 
with an assumed KATRIN prior for m^^ = 0.3eV. 



arguments as in the case where we used WMAP data only. A difference can be seen in the 
effect on n^. When using only WMAP data, a larger m^^ pulled Us to lower values, while 
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here Ug is shifted to shghtly larger values. This can be understood by the effect of M^, on 
the matter power spectrum, where a larger M,^ suppresses small scale fluctuations. This 
can be compensated by increasing n^. 

4-2. Constraining w 

Next we redo the analysis, including w as a free parameter. The resulting constraints on 
w can be seen in Figure 4 and Table 3. 




Figure 4. Marginalized parameter distributions of w when using data from WMAP 
only (left panel) and WMAP+LSS+SNla+HST+BBN+BAO+CMF(right panel). Black 
solid lines show the distributions for cosmological data only. This is compared to the 
resulting distributions when adding KATRIN data with m^^ = OeV (red dashed lines) 
and m^^ — 0.3eV (blue dotted lines). 



Cosmological data 


m^^ = free 


m^^ = OeV 


m^^ = 0.3eV 


WMAP 
WMAP++ 


1 q+0.5 
-'-•'J-0.6 

-1.06l°:J? 


-1.1 ±0.4 
-1.01 ±0.07 


-1.3 ±0.5 
-1.18 ±0.09 



Table 3. Limits on the equation of state of dark energy, w, with different 
priors on the neutrino mass. The first row are the results from using WMAP 
data only. The second row comes from using the full range of data sets 
WMAP+LSS+SNla+HST+BBN+BAO+CMF. 

When using WMAP data only, the constraints on w are relatively weak. One usually 
needs to include e.g. SNla or LSS data to get tight constraints on this parameter, as the 
main effect of w is to change the expansion history at late times. We see, however, that 
the lower limit on w improves from w > —2.5 to w > —1.8 (95% CL.) when including the 
KATRIN prior for m,,^ = OeV. This happens since increasing w will decrease the late-time 
expansion rate of the universe and thus shift the peaks in the CMB power spectrum to 
the left. This can be compensated by reducing M^, which will shift the peaks back to the 
right. 

When using the full range of cosmological data sets, the results become more 
interesting, as the constraints on w are tighter in this case. Here, the most interesting 
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effect occurs in the case of adding a KATRIN prior for rriy^ = 0.3eV, wliicfi makes w = — 1 
excluded by approximately 2a. That means that a positive KATRIN detection of ttIi,^ in 
this range, would be an indication of w < —1, which would be very interesting from a 
cosmological point of view. 

The reason for this dependence on M,^, can partially be explained by the shift of the 
peaks in the CMB power spectrum. But maybe more important is the fact that a smaller 
w will require a smaller Q\ to accommodate the accelerated expansion observed in the 
SNla data. This will in turn increase flm, which will allow for larger neutrino masses, as 
explained earlier. 

5. Goodness of fit 

Another interesting issue is the change in when adding the KATRIN priors. For each 
combination of cosmological model and data sets, we have calculated the change in by 

Ax^ = 2 In >C(no prior on M^) - 2 In >C(KATRIN prior on M^) , (5) 

where C is the likelihood of the best-fit sample in each combination of data sets on 
cosmological model. Thus Ax^ is a measure of how much worse the fit to the data becomes 
by introducing the KATRIN priors in the different cases. The resulting values of Ax^ are 
summarized in Table 4. 

We see that in the case of rrii,^ — OeV, the changes in x^ are small compared to 
the models with no prior on M,y. This is not surprising, as values of M^, ^ OeV fit the 
cosmological data very well. When using the full range of data sets in our 7 parameter 
model, the situation becomes slightly worse, giving Ax^ = —4.08. Introducing w as a 
free parameter in our model, the change in x^ by introducing the mj,^ = 0.3eV is not 
that severe anymore, giving Ax^ = —1.59. This can be understood by the well-known 
degeneracy between M^, and w. 



Cosmological data 


Model 


Ax\m,^ = OeV) 


Ax'(^.e = 0.3eV) 


-^par 


WMAP 


ACDM + 


-0.01 


-0.18 


7 


WMAP++ 


ACDM + 


-0.02 


-4.08 


7 


WMAP 


ACDM + M^ + w 


-0.13 


-0.01 


8 


WMAP++ 


ACDM + M^ + w 


-0.07 


-1.59 


8 



Table 4. Ax^ of the models with a KATRIN prior on M^, relative to the models 
with no prior on M^. WMAP -|— I- refers to the analysis with using all data sets 
WMAP-FLSS-l-SNla-FHST-l-BBN-FBAO-l-CMF as described in the text. iVpar refers to 
the number of free parameters in the models. 



6. Discussion and conclusions 

In this paper we have investigated whether constraints on m,^^ from the KATRIN 
experiment will affect our knowledge on cosmological parameters. This has been done 
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for two scenarios, one where rriy^ = OeV, and one where m^^ = 0.3eV. We have carried 
out the analysis both with a simple 7-parameter model with a cosmological constant, and 
extending the parameter space to include the equation of state for dark energy, w, as a 
free parameter. 

When using WMAP data only, we find that knowledge from the KATRIN experiment 
will contribute significantly to constrain a wide range of cosmological parameters, regardless 
of which of the m^^ scenarios we use. For instance will the significance of < 1 depend 
on what KATRIN tells us about m,y^. Other parameters that are sensitive to the value of 
rrii,^ are Q^rn, ^a, cts and Hq. 

Adding more cosmological data sets, both from SNla, galaxy catalogues and other 
priors, the situation changes a bit. In this case is strongly constrained from above by 
cosmology alone, such that an additional KATRIN prior in the case of rrii,^ — OeV has 
httle effect on our cosmological parameter constraints. However, if KATRIN measures a 
neutrino mass of rrii,^ — 0.3eV, there will be significant shifts in several of the parameter 
distributions. One should also note that several of the extra cosmological data sets added 
here may be affected by uncontrolled systematics (see [8]). Therefore, having cosmological 
constraints from WMAP+KATRIN without any additional cosmological data sets will be 
interesting regardless of the possibility to add other cosmological data sets to obtain similar 
results. 

In the case of w the most interesting result occurs in the scenario of a KATRIN 
defection of m,^^ = 0.3eV and using the full range of data sets. In this case, w < —1 is 
favored at a 2a level. It should also be mentioned that there are degeneracies between 
parameters from different cosmological inflation models and neutrino masses (see [33]). 
This means that a KATRIN prior on M^, will be important also for constraining inflationary 
models. 

To conclude, we find that the expected limits on rrii,^ from KATRIN, will be a useful 
input to constrain cosmological models, regardless of the value of rrii,^. If KATRIN detects 
a non-zero value of m^^ , this would be especially interesting. 
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